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The energy spectrum of ultra-high energy cosmic rays (UHECRs) is discussed on the basis of the Yakutsk array 
database analysis. In the region Eo = 10^'' to 3 X 10^^ eV the showers are detected under trigger-500, while 
at energies above 3 x 10^^ eV the whole acceptance area for trigger- 1000 is used in order to utilize all the data 
available in the region of Greisen-Zatsepin-Kuzmin (GZK) cutoff. Updated relations are adapted to estimate the 
primary energy having detected densities of charged particles at the distance 300/600 m from the shower axis. 



1. Introduction 

The primary energy estimation method is one 
of the clue features to deal with using extensive 
air shower (EAS) technique of ultra-high energy 
cosmic ray measurements. In this paper we dis- 
cuss the present status of algorithms used to as- 
sign the primary energy to showers detected with 
the Yakutsk array. Special attention is paid to 
the total flux measurement of the Cherenkov light 
on the ground which is the basis of the model- 
independent approach to the energy evaluation 
of EAS original particle. The uncertainty limits 
of the energy estimation method are analyzed. 

Resulting energy spectrum of cosmic rays in the 
region surrounding GZK energy Egzk ~ 4 x 10^^ 
eV has the form consistent with the prediction of 
Greisen pP and Zatsepin and Kuzmin |2]. One 
more EAS event at energy above 10^° eV has 
been detected with the Yakutsk array 18 Febru- 
ary 2004. An updated lateral distribution func- 
tion (LDF) of charged particles was used in data 
analysis which leads to increased Sgoo parameter 
in average. As the result we have 3 showers of 
energy greater than 10^" eV detected in Yakutsk 
to present day. 
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2. Data acquisition and analysis 

Actually, the Yakutsk array has two kinds of 
triggers to select the showers from the back- 
ground: produced by Cherenkov light detectors 
and by the scintillators. We will describe here 
the latter one; in the case of coincident signal 
(particle density p > 0.5 m~^) in two scintillators 
of each station within 2 microseconds the signal 
passes on to the central controller. Trigger-500 
is then produced in the case of coincident signal 
(in 40 /is) from three or more stations at ~ 500 
m spacing. Similarly, trigger-1000 is produced by 
~ 1 kilometer spacing stations. After 1992 when 
18 new stations were added, the array area is in- 
creased from 2.5fcm^ to 7.2km? where trigger-500 
operates. That is why we can deal with the spec- 
trum in the energy range from 3 x 10^^ eV to 
3 X 10^^ eV using the same trigger. 

In our previous analysis the lateral distribu- 
tion of charged particles was approximated by 
Greisen's formula "S*: 

f{r)^{-)-'\l + -f-'-\ (1) 
ro ro 

where b = 1.38 + 2. 16 cos 6* + 0.15 logio 5*600; 

here we use another approximation given ear- 
lier [ini: 

fir) - (-)-'-'(l + -)^-='-''(l + :-^)-^ (2) 

7-0 Tq zOOO 

where parameter c has the value given in Tabled 
and b is additionally adjusted; 
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Table 1 

Lateral distribution index 



logio£:o < 18.0 18.1 18.3 > 18.7 
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which has been shown to describe better the ex- 
perimental points at large core distances r > 1000 
m. As a result, the charged particle density at 600 
m from the core, 5*600, is increased by ^ 10% for 
showers above 2 x 10^^ eV with axes inside array 
perimeter, and by ^ 20% outside. 

3. Energy estimation 

The main distinctive feature of the Yakutsk ar- 
ray is the air Cherenkov light measurement. The 
total flux of the light emitted in atmosphere, Qtot, 
is used as the main estimator of the primary parti- 
cle energy. In order to derive the relation between 
Qtot and ionization loss of the shower electrons in 
atmosphere we have to use the classic formula by 
Tamm and Frank ,5, for the number of Cherenkov 
photons in the wavelength interval (Ai, A2) emit- 
ted when relativistic electron moves in air I cm 
with the speed v: 



a{x, E) = 2'Kal{ 



1 



)(1 



(3) 



where a = 1/137; n is refraction coefficient. 
Emission is possible when the particle velocity 
is greater than c/n. At the depth x in atmo- 
sphere the threshold energy of electrons is Ethr = 
rac? I ^j2{nQ — l)p{x)/ p^, where rip, po are air re- 
fraction coefficient and density at the sea level. In 
the first approximation with regard to rru? jE the 
number of photons emitted on the length \ g / cm?' 
is 

a(x,i?)« 4^(1-^)^(1-:^). (4) 
A2 Ai Po E^ 

The total light flux at the sea level is given as the 
integral of electron energy spectrum in the shower 
f{Ea,E,x) = -dNjdE: 

Qtot{Eo, Xq) — 

= - dxk{x,xo) dEf{Eo,E,x)a{x,E), (5) 

Jo J Ethr 



where k{x,xo) is attenuation coefficient of the 
light intensity in atmosphere. Integrating 
by parts and using the boundary conditions 
a{x, Ethr) — 0, Ne{EQ,EQ, x) = we take 

QtotiEo,XQ) = 
dxk{x,xo) / dENe{Eo,E,x)a{x,E). (6) 

J Ethr 

Due to slow variable Ne{Eo, E,x) at Ethr 
and a'{x,E) oc E^y^^/E^, we can factor 
Ne{Eo, Ethr, x) outside the integral over energy: 

Qtot{Eo,xo) « 47ra(-!- - -^) x 

X / dxk{x,XQ)Ne{Eo, Ethr,x). 

Po Jo 



(7) 



The energy spectrum of electrons near the shower 
maximum gives Ne{Eo, Ethr, x) = Ne{Eo,0,x) x 
77, where rj — 1 — ee*^ dte~^ /t is independent of 
X] e — 2.29Ethr/ P- The cascade curve of electrons 
can be approximated according to the measure- 
ments by HiRes [H]: 

Ne[Eo,0,x) ^ Nmaxexpi -[ — )), 



(7 X H~ '^iXjjiax 

where a = 0.272 ± 0.002. Attenuation of the 
Cherenkov light in atmosphere has the length 
A [xq — x) 4 . Substituting exp{{xo — a;)/ A) to 
the total flux integral and expanding near Xmax 
we have 

Qtot{Eo,xo) « 47ra(-!- - -^)— x 

A2 Ai Po 

xvHx^ax,xo){l + 5Q) / dxNe{EQ,0,x), (8) 
Jo 

\6Q\ < 5%. On the other hand, Ei = 
f3/to Jq" dxNe{Eo,0, x) is ionization integral. 
The final relation is 

Qtot{Eo, xo) « 47rQ(-!- - ^) x 
A2 Ai 



no ~ I .toEi 

X "TlfiyXmax, Xq)' 



with the accuracy ~ 5%. It is determined by 
Xmax and attenuation of light in atmosphere but 
is independent of hadronic interaction model as- 
sumptions. 
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Hereafter we analyze the shower parameters 
governing the energy fractions transferred to EAS 
components. 

3.1. Energy balance of EAS components 

The energy fractions of EAS primary parti- 
cle transferred to the shower components can 
be described on the basis of hadron transport 
equations j4l7| . If Ek, {k = N.ir.jiv.e^) is the 
energy transferred to nucleons, charged pions, 
muons+ncutrinos, electrons+photons, then we 
are going to demonstrate in this subsection, with 
simple arguments, that a few cascade parame- 
ters determine the ratios between E^. - the energy 
balance in the shower. For instance, the trans- 
port equation for charged pions density 7r(a;, E) 
at depth x: 

dTT{x,E) 1 

+ ^ / " Ax,U)w^^iE,U)dU + 

+^ f " N{x,U)w^N{E,U)dU, (10) 

where interaction mean free paths A^r, Ajv are as- 
sumed constant; 'WT^T^{E,U),WTtN{E,U) are the 
spectra of charged pions produced in pion-air and 
nucleon-air interactions, can be transformed (in- 
tegrating J EdE) to: 



E^ B^tt{x, E >0) 2E^ 2KnE_ 



N 



(11) 



dE^ 

dx Xtt X ' SAtt SAat 

where E„{x) = J^^" tt{x, E)EdE; tt{x,E > 0) = 

J^" 7r(x, E)dE; Kn is nucleon inelasticity sup- 
posed to be constant; En = Eoexp{—KNx/XN)- 
In the energy range E » B^^ the only pa- 
rameters to define the quantity Et^ are Km/Xn 
and Att. It means that the energy transferred to 
charged pions is independent of the spectra of pi- 
ons produced in nuclear interactions. Hence, in 
the general case, we can use the simple i5-model 
with Wik{E, U) = nsS{E — U /us), where Ug is the 
multiplicity of secondaries, to balance the com- 
ponents energy in a shower. Because of the net 
value of Efj^+^/Eo < 0.1, the uncertainty due to 
simplified model should be of the second order of 
magnitude. 



To summarize, the main model parameters to 
govern the energy balance in the shower are aver- 
age inelasticity coefficients, mean free paths, mul- 
tiplicity of secondaries and the fragmentation rate 
of primary nucleus. Other model characteristics 
such as 'the form of rapidity distribution of con- 
stituent quarks' are of less influence. The resul- 
tant energies in the case of constant A^, Aat, Kn 
and = are shown in Fig. ^ together with 
(5-modcl results. An asymptotic (x = oo) estima- 
tion of Eej with CORSIKA(+QGSjet) program 
at Eq = 10^^ eV (open circle) is shown here as 
well. 




400 600 800 

X, g/cm2 



Figure 1. The energy carried by the cascade nu- 
cleons (N) , charged pions (tt) , muons and neutri- 
nos {fi+iy), electrons and photons (e-l-7). Dashed 
curves show analytic expressions with A^ , — 
const, Bt:: = 0: En = Eoexp{—KNx/ Xn)', 
Et; = 2/3i?o(l — exp{—KNx/XN))exp{—x/XTr/'i)', 
Ee-y — Eq — En — Ej^. Solid curves are S- 
model results: cross sections are supposed ris- 
ing oc O.OS^n^;; B^ = 120 GeV; Kn = 0.5; 
n, oc E^/^; Eq = 10^* eV. 



3.2. Experimental evaluation of the energy 
dispensed to EAS components 

Energy fractions of the main EAS components 
can be estimated using the Yakutsk array data. 
Ionization loss of electrons is measured here de- 
tecting the total flux of the Cherenkov light on 
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Table 2 

The primary energy portions gone with EAS com- 
ponents. Eq = 10^^ eV. 9 = 0°. Notations: 
is ionization loss of electrons in the atmosphere; 
Eg is ionization loss of electrons in the ground; 
Efi+u is energy transferred to muons and neutri- 
nos; Eh is energy carried by the nuclear active 



component. Eq ~ 


Ei + Eg + £"^4 


+ Eh. 


Energy deposit 


The portion 


Experimental 


channel 


of energy, % 


uncertainty, % 


E^ 


78 


30 


Eg 


9 


30 




9 


16 


Eh 


4 


20 



the ground. A relation between these values is 
given by formula lO with parameters relevant to 
the Yakutsk array, taking into account detector 
calibration and atmospheric transparency factors 

^ ^ 2.18 X IQ^Qtot 

" 0.37-1-1.1 X 10-3X„aa;' 

Other portions of the energy carried out by 
electromagnetic and muonic components beyond 
the sea level is evaluated via the total number of 
electrons: 

Eg — eoNeXe/to, 

where eo , are the critical energy and radiation 
length of electrons in air; attenuation length Ae is 
derived from zenith angle dependence of N^, \^ ; 
and muons measured on the ground: 

E^ = N^{E>1 GeVW;, 

where the average energy of muons E'p is taken 
from the MSU array data [T^ . 

Residuary energy fractions transferred to neu- 
trinos -Ey, nucleons Eh etc., unmeasurable with 
this array, are estimated using model calculations. 
The resulting apportioning of the primary energy 
10^8 eV is given in Tabled 

The energy fraction carried by electromagnetic 
component appears to be the basic contribution 
to the total energy of the shower, and its energy 



dependence (measured with Cherenkov light de- 
tectors -t- scintillators of the Yakutsk array) is 
illustrated in Fig. |21 together with (5-model es- 
timation under superposition assumption for the 
primary proton and iron nucleus. 
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Figure 2. The electromagnetic component en- 
ergy estimation from the Yakutsk array data: au- 
tonomous sub-array (rhombuses), the main array 
(squares); P, Fe curves are (5-model results. 



Because the air Cherenkov light total flux, and 
electron and muon number of the shower are ex- 
perimental values measured on the ground, only 
about 10% of the primary energy Eq = 10^® eV 
is calculated using the model assumptions. So we 
consider the energy estimation algorithm in use 
in the Yakutsk group to be model-independent in 
the first approximation. 

Moonless nights when air Cherenkov light mea- 
surements are possible constitute ~ 10% of the 
observation period. In order to evaluate the pri- 
mary energy of the bulk of showers, the linear 
correlation 

logio ^600 = l.OSlogio Q400 - 0.79; = 0.997 

is used between the charged particle density at 
600 m from the shower core and the light intensity 
at 400 m from the core, Q400 (Fig- El which, in 
turn, is related to the total flux of the Cherenkov 
light on the ground 
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Figure 3. Linear correlation between charged par- 
ticle density at 600 m from the core (iSeoo) s^nd air 
Cherenkov light intensity at 400 m {Q400) mea- 
sured in the same showers. 9 = lO". 



Experimental uncertainties in EAS component 
energies estimated using the Yakutsk array data 
are summarized in Table |21 The main contri- 
bution arise from 6Ei which is formed by un- 
certainties in atmospheric transparency (15%), 
detector calibration (21%) and the total light 
flux measurement (15%). Errors in estimation 
of iVe, Ae, Nfj_ + Ni, determine the next two items 
(for ionization in the ground and SE^j^,^). Resul- 
tant energy estimation uncertainty is the sum of 
all errors weighed with the second column of the 
Tabled SE^ - 30%. Extra 20% is added due to 
Sem — Q400 conversion uncertainty. 

The conversion relation from 6*300 / -S'eoo to pri- 
mary energy was re-examined in 2003 11 . For 
vertical showers {x = lQ2Qg/cm?) we have: 

Ea - (6.5 ± 1.6) X 10^6^300 (0°)°■^^±°■°^ 

Eo - (4.6 ± 1.2) X 10^^600 (0'')°-''±° °'. (12) 

Uncertainty of the conversion factor is mainly due 
to absolute calibration error of the Cherenkov 
light detector and optical transparency estima- 
tion error of the atmosphere. These errors affect 
the intensity of CR flux but not the shape of the 
spectrum. 

The observed densities S'3oo/5'6oo s-t zenith an- 
gle, 6*, are connected to the vertical direction val- 
ues along attenuation curve |12j . In order to mea- 



sure the attenuation length of these densities for 
fixed primary energy, we have used two diflerent 
methods - well-known equi-intensity cut method, 
and fixing the Cherenkov light intensity at 400 
m from the core as the equivalent of the pri- 
mary energy, taking into account the light absorp- 
tion in the atmosphere. In Fig. 0] the results are 
given. Experimental points are consistent with 
each other for the two methods used and can be 
described by the sum of two components - soft 
(electrons, attenuation length Ae — 200 g/cm?) 
and hard (muons, A^ = 1000 y/cm^): 

Szoo{e) = ^30o(0°)((l - /330o)ea;p((a;o - x)/K) + 

+l3'iOQexp{{xo - x)l\^j,)), (13) 

where f3zm is the hard component fraction. At- 
tenuation curve for 6*600 is the same but /? is dif- 
ferent: 

/?3oo = (0.563 ±0.032)53oo(0°)-°■l«^"■°^ 
/?600 = (0.62 ± 0.006)56oo(0°)-° °'^±° °^ (14) 

Previously we used to apply another attenua- 
tion curve in order to estimate the primary en- 
ergy 13] which leads to divergence in energies of 
inclined showers rising with zenith angle. 




Figure 4. 6*300 as a function oix = 1020/ cosO for 
diflerent CR intensities. Open symbols [Si] are 
equi-intensity method results, filled ones {Qi) are 
derived fixing (5400- 
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Four showers in the energy range Eq > 10^" cV 
are given in Table |31 (one event is added slightly 
below the threshold because the energy estima- 
tion error is larger than tiny difference). 



Table 3 

Largest EAS events detected with the Yakutsk 
array 



Date 


0" 


\ogEo 


5Eo,% 






18.02.04 


47.7 


20.16 


42 


16.3 


140.2 


07.05.89 


58.7 


20.14 


46 


2.7 


161.6 


21.12.77 


46.0 


20.01 


40 


50.0 


220.6 


15.02.78 


9.6 


19.99 


32 


15.5 


102.0 



20 
18 
16 
14 
12 
10 
8 

6 ^ 

4 
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■ 



17 



18 19 

log(E, eV) 
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4. The energy spectrum measured with 
the Yakutsk array 

To derive SeffTQ. - the product of array ac- 
ceptance area, observation periods duration and 
solid angle, we have used = tt; the real duration 
of the array duty time periods as T; and have cal- 
culated the area ^e// (5*300/600 1 ^) limited within 
the array perimeter. The array configuration has 
been changed several times since 1974; the ac- 
tual perimeter was used in particular time pe- 
riod. For the highest- energy showers Ei^ > Eqzk 
the extended perimeter was applied. In this case 
the array area is increased ^1.4 times while the 
shower parameters reconstruction is still possible 
with acceptable error. Energy dependence of the 
array area calculated for the vertical showers in 
the case of strict perimeter and extended one is 
shown in Fig. [S] and Fig. correspondingly. The 
same procedure is applicable at each zenith angle 
for inclined EAS because S{0) ^ S{0°) sec 9. 

In Fig. resulting integral energy spectrum of 
cosmic rays is given composed of showers above 
3 X lO^"^ eV, with zenith angles 9 < 60°. We 
confirm with the present dataset analysis the ex- 
istence of an 'ankle' feature in the shape of spec- 
trum around 10^^ eV, i.e. the flattening in the 
spectral index revealed earlier by all the groups 
working in this energy range T3,. At the high- 
est energy end of the spectrum we still have a 



Figure 5. Array area for vertical showers vs Eq. 
Circles are for the period before 1990; rectangles 
are for after the array was contracted. Two ap- 
proximation functions of lateral distribution of 
charged particles are used: formula ^ (filled 
signs) and formula Q (open signs). 



deficit in agreement with the expected flux cutoff 
due to interactions of CRs with relic microwave 
background. As an example of the model cal- 
culation results, the curve in Fig. [7| shows UHE 
proton spectrum |18| in the model with uniform 
distribution of extragalactic sources with Ei ^2 = 
6.2 X 10^^ eV. 

Comparison with other spectra measured with 
AGASA and HiRes arrays is shown in Fig. |S1 
Differential intensities in energy bins measured 
in Yakutsk and statistical errors of the data are 
given in Tables 01 El for usability of the spec- 
trum. In the energy range below 10^^ intensi- 
ties are highest of the Yakutsk array data. Most 
likely it is due to different energy estimation tech- 
nique used in three cases - resultant energies di- 
verge. In AGASA case the model dependent re- 
lation between Sqqq and Eq is used, HiRes team 
measures the fluorescence light flux which is con- 
nected to ionization loss of electrons in atmo- 
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Figure 6. The same areas but for extended 
perimeter. Notations are the same. 



sphere. Our relations 11211 for vertical showers 
give energy higher than model calculations by 
(30-40)% at iJo ~ 5 X lO" eV and by (15-20)% 
at £^0 ~ 5 X 10^9 eV. 

At the highest energies our intensities became 
closer to AGASA data, but GZK cutoff feature is 
different. The slope of the spectrum above 10^" 
eV seems rising in the case of the Yakutsk array 
data as well as for HiRes spectrum; on the con- 
trary, AGASA data exhibit no cutoff. But the 
probability to detect three or less events above 
10^*^ eV in Yakutsk is 15%, while 6 events are ex- 
pected according to AGASA's intensity. There- 
fore, it seems that one has to wait until Auger/TA 
data in order to decide whether there is the GZK 
cutoff or not. 

In addition to the charged particle detection 
on the ground, we have another technique at the 
Yakutsk array - the air Cherenkov light measure- 
ment, which can be used to draw out the CR 
spectrum in independent way These two 

spectra are shown in Fig. |3 covering the wide en- 
ergy region from 10^^ to 10^° eV. The Cherenkov 
detector sub-array is covering only about 3 km^ 



Table 4 
Cosmic ray 
within array 



flux measured 
perimeter. 



under trigger-500. 



Eo 


J(E) X E'^ 




oJ 


5.89E+19 


4.76E+24 


2.75E- 


1-24 


3.89E-hl9 


1.15E+25 


2.86E- 


1-24 


2.57E-hl9 


4.90E-f24 


1.19E- 


1-24 


1.62E+19 


4.20E-h24 


7.09E- 


1-23 


9.55E-hl8 


4.15E-I-24 


4.15E- 


1-23 


7.94E-hl8 


4.56E-I-24 


5.29E- 


1-23 


6.31E-)-18 


4.41E-h24 


4.14E- 


1-23 


4.90E-hl8 


4.69E-f24 


3.42E- 


1-23 


3.98E-hl8 


4.81E-h24 


2.79E- 


1-23 


3.16E+18 


4.89E-f24 


2.30E- 


1-23 


2.51E-hl8 


5.68E+24 


1.99E- 


1-23 


2.00E-hl8 


6.08E-h24 


1.70E- 


1-23 


1.58E+18 


5.68E-f24 


1.36E- 


1-23 


1.26E-hl8 


5.75E+24 


1.09E- 


f-23 


1.00E-hl8 


6.17E-h24 


1.05E- 


1-23 


7.94E-hl7 


6.84E-I-24 


1.03E- 


1-23 


6.31E-hl7 


7.05E-I-24 


1.06E- 


f-23 


5.13E-hl7 


7.74E+24 


1.39E- 


1-23 


4.27E-)-17 


7.64E-I-24 


2.75E- 


1-23 



Table 5 

Cosmic ray flux measured under trigger- 1000, 
within array perimeter. 



Eo 


J{E) X 




6J 


4.07E+19 


5.75EH 


h24 


1.29E- 


1-24 


2.88E-hl9 


5.08EH 


1-24 


8.59E- 


}-23 


2.09E+19 


4.23EH 


h24 


5.58E- 


^23 


1.48E-fl9 


4.70EH 


h24 


4.23E- 


F23 


1.02E-)-19 


3.47EH 


h24 


2.53E- 


f-23 



Table 6 

Cosmic ray flux measured under trigger- 1000, 
within extended array perimeter. 



Eo 


J{E) X E"^ 


SJ 


1.41E+20 


4.5GEH 


h24 


3.18E+24 


l.OOE+20 


2.3GEH 


h24 


1.628E+24 


6.43E-hl9 


7.19EH 


h24 


1.75E-I-24 


4.81E+19 


6.G4EH 


h24 


1.23E-h24 
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Figure 7. Integral energy spectrum of CR de- 
tected with the Yakutsk array. Data selected 
within array perimeter: open circles (trigger- 
500); ominus (trigger- 1000); filled circles (trigger- 
1000, extended perimeter). 



Figure 8. Differential spectrum of UHECR, mul- 
tiplied by Eg. For the Yakutsk array data nota- 
tions are the same as in Fig. Open triangles 
show the Akeno ^ result, filled ones - AG AS A 
|15|. Open rhombuses are for HiRes-I data and 
filled ones for HiRes-II [T7|. 



area, so it is convenient at energies E < Egzk, 
while scintillator stations can measure the spec- 
trum above this limit. At the low energy thresh- 
old of the main array ~ 10^^ eV we have a dis- 
crepancy in intensity measured using two meth- 
ods. This may be caused by S^qq uncertainty in- 
creasing in the vicinity of the threshold. 

Two features of the measured spectrum - a 
'knee' (at 3 x lO^^ eV) and 'dip' around 

E 10^^ eV can be explained in the model 
of anomalous diffusion of galactic cosmic rays in 
fractal interstellar medium ■ In this model 
the mass composition is predicted in the energy 
range 10^'' — I0^° eV; EAS maximum depth vs. 
energy is simulated by CORSIKA v. 6.0 code us- 
ing QGSjet model assumptions. Resulting spec- 
tra of the primary nuclei groups are shown by 
curves in Fig. 13 together with the total spectrum 
marked 'AH'. 

5. Conclusions 

We have measured the energy spectrum of cos- 
mic rays below and surrounding the GZK cutoff 
region. Air Cherenkov light measurements pro- 



vide us the model-independent approach to the 
primary energy estimation using the total light 
flux connected to the ionization loss of electrons 
in atmosphere. Other appreciable energy frac- 
tions are evaluated on the basis of the total num- 
ber measurement of electrons and muons on the 
ground. An updated function to approximate the 
lateral distribution of charged particles leads to 
increased average density Seoo- As a result, we 
have three showers detected with energies above 
10^" eV including a new event detected 18.02.04. 
Two energy spectra measured with the Yakutsk 
array cover the range 10^^ to ~ 10^" eV and 
exhibit the 'knee' and 'dip' features in agree- 
ment with other arrays. The spectrum tail above 
Egzk is consistent with the GZK cutoff contrary 
to AG AS A result, although the statistical signif- 
icance of this conclusion is insufhcient due to a 
few showers detected above 10^° eV. 
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